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[ We analyse the motion of solitary type optical pulse in a fused optical fibre coupler 
when inter channel raman scattering (IRS) is taken into account and same type of 
dispersion management is used for both the core. The paper is divided into two 
sections . In the first half we utilise a moment method to study the variations of 
the different parameters of the Gaussian pulse, such as chirp, frequency, energy etc. 
while in the second half assuming a general form of the stationary pulse we discuss 
the switching phenomena due to the coupling. It is shown that lauching a pulse in 
one arm of the coupler can generate pulses in both the arms, but the patterns differs 
in a significant manner when the strength of IRS gets changed. ] 
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1 Introduction 



Nonlinear directional couplers are now a days the building blocks of all optical pro- 
cessing. They are now subject of intense theoretical and experimental research[l]. 
The use of soliton like pulse in conjuction with these couplers has brought a revolu- 
tion in the field of nonlinear optics [2]. Most directional couplers are usually twin core 
coupler with a parmanent coupling between the cores [3]. On the other hand inter- 
channel Raman scattering is an important phenomenon which plays a decisive role 
in the formation and propagation of a pulse in a fibre [4]. So, in this communication 
we have considered the mutual effect of intrapulse raman scattering (IRS) and the 
coupling on the propagation characteristics of the pulse in a fused coupler. These 
types of couplers are now a days very important as they can be used for switching 
operation[5] . Here we have observed that a very weak pulse launched in one arm of 
the coupler can grow considerably due to the coupling[6]. A similar event takes place 
in case of conversion, compression and splitting of pulses. 

The paper is organised as follows. In the first section we derive the equations 
governing the evolution of the pulse parameters such as energy, width, centre, chirp 
etc. by following a moment method[6a]. Because of the presence of IRS, a variational 
procedure depending on a Lagrangian cannot be used[7]. These are then numeri- 
cally solved and these parameters variations are studied with respect to the coupling 
strength and the magnitude of IRS. 

In the second half of the present communication we have considered a general 
form of the stationary pulse and studied the effect of the coupling in the generation 
of pulse in both the arms when there is only one pulse initially in one arm. It is 
observed that the presence of IRS strongly effects such switching pattern. 

2 Gaussian pulse parameter evolution: 

The equation describing the propagation of the nonlinear waves in a fused coupler is 
written as, 

du i3 2 d 2 u ( . l9 _ d , ,,\ . , 
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where u(z,t) and v(z,t) represent the wave forms in the two arms of the coupler, 
^{z) is the dispersion management whose explicit forms are given in equations (20) 
and (21), T R is the coefficient of interchannel Raman scattering and K is the coupling 
constant. This coupled set is difficult to solve but can be studied numerically. One 
of the most popular method is to assume an initial profile of the incomming wave 
form (outgoing) and to study the variations of the physical parameters of the pulse 
propagating along the fibre. Unfortunately, the presence of the IRS term does not 
allow a suitable Lagrangian for the system. So the standard variational technique can 
not be implemented and we take recourse to the moment method. The basic idea is 
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to treat the optical pulse like a particle whose energy (E), position (T), and frequency 
shift (ft) (from the original carrier frequency) are defined as 
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The root mean square (RMS) width of the pulse is, 
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The chirp of the pulse is given as, 
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where ip stands for either u(z,t) or v(z,t). The initial profile of the pulse u,v are 
assumed to be Gaussian-type, 
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Differentiating (3) to (7) and using (1), (2) and (8) we get the following system of 
nonlinear coupled evolution equations , 
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This set of ten ordinary differential equations is then integrated by the Runge-Kutta- 
Felberg algorithm. In the above equations, 
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The above ordinary differential equations can now be used to simulate two different 
physical phenomena. One is the coupling between the two modes of propagation on 
the two cores and the other is the effect of Raman interchannel scattering and not 
the least their mutual influence on each other. We also show the initial and final form 
of \u\ 2 and \v \ 2 to study the effect of switching in the coupler. Before proceding to 
the actual numerical simulation, we comment on the type of dispersion management 
used, 
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Type (a) 

with dispersion map, 
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and map strength, 
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Type (b) 
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The type (a) was used by many workers[8] with an eye to the actual construction of the 
anomalous and normal section of a optical fibre, with specified gap of amplification. 
On the other hand the second type of dispersion management was used by Serkin 
etal[9]. for explicit analytic construction of soliton through inverse scattering to be 
used as a model of optical pulse. 

The results of our numerical simulation are clearly displayed in figures. To start 
with consider figure(l) where the variation of E(energy), Q(RIFS), r(width) and 
C(chirp) are plotted with respect to distance travelled by the pulse. In each of the 
figure we have shown the effect of variation of the coupling constant K keeping all 
other parameters fixed. The type of dispersion management used is given in equation 
(21), with S = —0.9. On the other hand in figure (2) the same type of dispersion 
management with different value of 5 (=0.9) is shown. In both the cases the value of 
T R is the same. On the other hand due to the usual type of dispersion map given 
in equation (20) , the variation of E, C, RIFS, r are exhibited in figures (3), (4), (5) 
for different values of K. One may note that the variation of energy (E) in figure (5) 
is shown that while it decreases in one core, in the second one it grows. The same 
phenomenon is clearly visible in figure (6) where the growth in one and decay in the 
other are exhibited. So this kind of phenomenon is actually a switching action for 
which such couplers are constructed. The final shape of output pulse in both the cases 
are exhibited in figures(6) and (7). The corresponding parameter values are given in 
the figure caption. It may be noted that even if we keep the amplitude of the input 
pulse in one of the channel very small, still there is a significant display of the pulses 
in both the cores. One may observe that it is not possible to set one of pulses exactly 
equal to zero due to some mathematical inconsistencies It may be added that all these 
analysis can easily be repeated if the initial profile of the pulse is assumed to be of 
such type, rather than Gaussian are. 
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3 Switching Action: 



The switching action can be further studied if we neglect the dispersion management 
and consider the simplified model which exhibits translational invariance. Assuming 
propagating pulse of the form 

u = f(z — out) and v = g(z — out) 

we get 

d£ + fde + = 11 1 ' u 1 + u u ~di^ u ' v ( ^ 
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Here ^ and 7 all are constants. Writing u = u\ + m 2 and i> = u 1 + iv 2 and separating 
real and imaginary parts we get four ordinary differential equations of second order. 
Assuming that at the launch condition v = 0, and u = uq , we have integrated the 
system by Runge-Kutta, ensuring the boundary condition needed for a solitary pulse. 
The shape of the pulse as it imerges is shown in the figure (8) for various coupling 
constan K. Note that here the input pulse is kept zero in one arm, but in the output 
side, it can be seen in both the arms. . 

4 Conclusion: 

Here in this note we have analyzed a model of a fused coupler, where both the 
dispersion management and IRS are assumed to be present. It is very clear that the 
type of dispersion management significantly effects the variation of the parameters, 
on the other hand the IRS causes some irregularities in their variation. 
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7 Figure Caption 



Figure-l. Variation of parameters, E, f2, r and C of u and v with respect to z taking 
5 = — 0.9where dotted(..) line stands for k=-3, dash-dotted(-.) line for k=3 and 
solid(-) line for k=0.9 and T R = 0.696, 7 = 0.1994 

Figure-2. Variation of parameters, E, Q, r and C of u and v with respect to z taking 
5 = 0.9where dotted(..) line stands for k=-3, dash-dotted(-.) line for k=3 and solid(-) 
line for k=0.9 and T R = 0.696, 7 = 0.1994 

Figure-3. Variation of E and C of u and v with respect to propagation distance when 
T R = 0.5, 7 = 0.1994 and k=1.0 

Figure-4. Variation of E and C of u and v with respect to propagation distance when 
T R = 1.5, 7 = 0.1994 and k=2.0 

Figure-5. Variation of flandr of u and v with respect to propagation distance when 
T R = 0.696, 7 = 0.1994 and k=0.09 

Figure-6. Output pulses with respect to propagation distance when T R = 0.696,7 — 
0.1994 and k=0.09 

Figure-7. Output pulses with respect to propagation distance when T R = 0.696,7 = 
0.1994an<iA; = 0.09 with different initial profile. 

Figure-8a. Output pulses with respect to propagation distance with 7 = 0.02, (5 2 — 
-0.02,u = 0.9, T R = 1.6 

Figure-8b. Output pulses with respect to propagation distance with 7 = 0.02, /3 2 = 
-0.02, w = 0.9, T R = 0.5 

Figure-8a. Output pulses with respect to propagation distance with 7 = 0.02, fi 2 — 
-0.02, w = 0.9, T R = 3.0 
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